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Abstract 

High  energy-density  plasma  (HEDP)  is  created  using  a 
plasma- filled  rod  pinch  diode  (PFRP)  [1]  on  the  Gamble 
II  generator.  HEDP  is  produced  by  a  1-2  MeV,  500  kA 
electron  beam  that  deposits  energy  at  the  end  of  a 
tungsten  rod.  When  the  rod  is  tapered  to  a  point,  the 
energy  density  is  estimated  to  reach  2.4  MJ/cm^  in  about 
10  ns,  before  rapid  expansion  decreases  the  energy 
density. 

Recent  measurements  show  that  the  temperature  peaks 
at  about  30  eV  at  the  time  of  maximum  energy  density, 
and  that  the  time-averaged  ionization  is  about  +17,  similar 
to  MHD  model  predictions  [2].  The  plasma  mass 
distribution  is  inferred  from  x-ray  distribution 
measurements.  The  time-dependent  mass  distribution  is 
used  to  calculate  the  mass  density  and  plasma  pressure. 
These  measurements  indicate  the  HEDP  is  strongly 
coupled  (E  =  35)  with  peak  pressure  (16  Mb)  about  7 
times  greater  than  ngkT. 


I.  INTRODUCTION 

High  Energy  Density  Plasmas  (HEDP),  generally 
defined  as  plasmas  with  energy  densities  greater  than  0. 1 
MJ/cm^  or  pressures  greater  than  1  Mb,  are  subjects  of 
current  research  [3].  One  challenge  is  to  produce  and 
diagnose  HEDP  in  the  laboratory.  This  paper  reports 
measurements  of  HEDP  parameters  created  in  tungsten 
using  a  Plasma-Filled  Rod-Pinch  diode  (PFRP)  [1]. 

The  PFRP  is  connected  to  the  Gamble  II  pulsed  power 
generator  (2  MV,  0.7  MA,  60  ns  with  a  matched  3 -Ohm 
load)  as  shown  in  Fig.  la.  The  positive- voltage  electrode 
is  terminated  in  vacuum  with  a  1  mm  diameter  tungsten 
rod.  The  end  of  the  rod  is  tapered  to  a  more-or-less  sharp 
point.  The  rod  extends  25  mm  through  a  hole  in  a  flat, 
grounded  cathode  plate.  Six  cable  guns  inject  plasma 
Starting  about  1  ps  prior  to  the  start  of  the  generator 
current.  Figure  la  is  a  notional  picture  of  the  plasma  after 
the  current  has  increased  to  about  500  kA  in  50  ns,  and 
after  being  displaced  in  the  axial  and  radial  directions  by 
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JxB  forces.  When  the  plasma  is  swept  away,  a  high- 
energy  (-  MeV)  electron  beam  connects  the  plasma  to  the 
rod  tip.  The  electron  beam  heats  the  tungsten,  creating 
HEDP,  and  produces  bremsstrahlung. 

Typical  electrical  and  x-ray  waveforms  for  the  PFRP  on 
Gamble  II  are  shown  in  Fig.  lb.  The  current  waveform 
(blue  line)  increases  to  over  400  kA  before  the  x-ray 
signal  (green  line)  starts.  The  voltage  (red  line)  between 
the  cathode  disk  and  the  rod  is  calculated  using  an 
upstream  voltage  monitor,  and  subtracting  LdEdt  where  L 
is  the  inductance  between  the  voltage  probe  and  the 
cathode  plane.  The  voltage  exceeds  1  MV  during  the  rise 
time  (about  10  ns,  yellow  shaded  region)  of  the  x-ray 
signal  and  the  coupled  energy  (magenta  line)  is  about  10 
kJ  in  this  time  interval.  The  zero  time  is  chosen  to  be  the 
Start  of  the  x-ray  signal  for  synchronizing  data  from 
different  shots. 

Previous  modeling  [2]  assumed  the  tungsten  plasma  had 
a  time-dependent  Gaussian  radial  profile  and  a  fixed 
length  of  3.5  mm,  consistent  with  time-integrated  x-ray 
images.  A  self-similar  hydrodynamic  expansion  model 
based  on  Sesame-tabulated  equation  of  state  data  indicates 
the  maximum  energy  density  is  about  2.4  MJ/cw?  at  10 
ns.  At  this  time,  the  predicted  tungsten  ionization  is  +10, 
T  =  25  eV  and  the  mass  density  on  axis  is  close  to  solid 
(19  g/cm^). 


II.  HEDP  DIAGNOSTICS 

The  diagnostics  used  to  measure  the  HEDP  plasma 
parameters  are  described  in  Fig.  2.  Filtered  x-ray  diodes 
(XRDs,  Fig.  2a)  are  used  to  measure  the  plasma 
temperature.  L-line  spectroscopy  (Fig.  2b)  provides  data 
that  determines  the  time-integrated  ionization.  An  array 
of  20  pin  diodes  (Fig.  2c)  records  x-ray  signals  behind  a 
thick  “rolled  edge”  to  measure  the  x-ray  distribution  as  a 
function  of  time.  The  x-ray  distribution  measurements 
determine  the  mass  density  distribution,  used  to  calculate 
the  plasma  pressure  (independent  of  temperature  and 
ionization). 


978-1 -4673-51 68-3/1 3/$31. 00  ©2013  IEEE 


316 


Report  Documentation  Page 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 


1.  REPORT  DATE 

JUN  2013 


2.  REPORT  TYPE 

N/A 


3.  DATES  COVERED 


4.  TITLE  AND  SUBTITLE 

High  Energy-Density  Plasma  Dynamics  In  Plasma-Filled  Rod-Pinch 
Diodes 

6.  AUTHOR(S) 


7.  PEREORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Plasma  Physics  Division,  Naval  Research  Laboratory  Washington,  DC, 
USA 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


8.  PEREORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES)  10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADM002371.  2013  IEEE  Pulsed  Power  Conference,  Digest  of  Technical  Papers  1976-2013,  and 
Abstracts  of  the  2013  IEEE  International  Conference  on  Plasma  Science.  IEEE  International  Pulsed  Power 
Conference  (19th).  Held  in  San  Francisco,  CA  on  16-21  June  2013.,  The  original  document  contains  color 
images. 

14.  ABSTRACT 

High  energy-density  plasma  (HEDP)  is  created  using  a  plasma-fllled  rod  pinch  diode  (PFRP)  [1]  on  the 
Gamble  II  generator.  HEDP  is  produced  by  a  1-2  MeV,  500  kA  electron  beam  that  deposits  energy  at  the 
end  of  a  tungsten  rod.  When  the  rod  is  tapered  to  a  point,  the  energy  density  is  estimated  to  reach  2.4 
MJ/cm3  in  about  10  ns,  before  rapid  expansion  decreases  the  energy  density.  Recent  measurements  show 
that  the  temperature  peaks  at  about  30  eV  at  the  time  of  maximum  energy  density,  and  that  the 
time-averaged  ionization  is  about  +17,  similar  to  MHD  model  predictions  [2].  The  plasma  mass 
distribution  is  inferred  from  x-ray  distribution  measurements.  The  time-dependent  mass  distribution  is 
used  to  calculate  the  mass  density  and  plasma  pressure.  These  measurements  indicate  the  HEDP  is  strongly 
coupled  (A  =  35)  with  peak  pressure  (16  Mb)  about  7  times  greater  than  nekT. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIEICATION  OE: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OE 

18.  NUMBER 

ABSTRACT 

OE  PAGES 

SAR 

6 

19a.  NAME  OE 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Time  (ns) 


Figure  1.  (a)  PFRP  setup:  MeV  electrons  (green  arrows)  strike  the  tip  of  the  tapered  tungsten  rod,  creating  an 
intense  x-ray  source  and  high  energy  density  plasma,  (b)  Electrical  and  x-ray  signals. 


A.  XRD  measurements 

The  XRD  technique  in  Fig.  2a  is  a  standard  method  for 
measuring  temperatures  in  the  10-1000  eV  range  [4]. 
Photons  from  the  PFRP  pass  through  a  filter  before 
interacting  with  an  aluminum-foil  photocathode.  Two 
XRDs  with  different  filters  (1.8  pm  polycarbonate  and  0.8 
pm  aluminum)  have  appropriate  responses  to  determine 
the  plasma  temperature  from  the  signal  ratio.  The 
polycarbonate-filtered  XRD  response  increases 
monotonically  while  the  aluminum- filtered  XRD  response 
decreases  monitonically  with  temperature  in  the  10-50  eV 
range.  Grids  (dotted  lines  in  Figl  2a)  reduce  the  photon 
flux  and  apertures  reduce  scattered  light  signals.  B  fields 
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Figure  2.  (a)  XRD  (temperature  diagnostic),  (b)  L-line 
spectrometer  (ionization  diagnostic),  (c)  PIN  diode  array 
(x-ray  distribution  diagnostic). 


from  permanent  magnets  deflect  charged  particles.  The 
distance  from  the  PFRP  to  the  photocathode  is  adjusted  to 
obtain  a  desired  signal  level.  On  Gamble  II,  the  distances 
were  1.3  m  for  the  polycarbonate- filtered  XRD  and  1.9  m 
for  the  aluminum- filtered  XRD. 

The  temperature  derived  from  the  XRD  signal  ratio  is 
shown  in  Fig.  3a  for  a  typical  PFRP  shot.  The  XRD 
signals  are  very  small  during  the  low  impedance  phase  of 
the  PFRP.  The  XRD  signal  begins  slightly  earlier  than 
the  hard  (bremsstrahlung)  x-ray  signal.  The  temperature 
reaches  its  peak  value,  30  eV,  during  the  rise  of  the  x-ray 
signal. 

B.  L-line  spectroscopic  measurements 

A  high  sensitivity  x-ray  spectrometer  [5]  records  time- 
integrated  line  spectra  from  the  PFRP  plasma.  The 
diagnostic  and  data  analysis  are  described  in  [6].  A 
typical  spectrum  recorded  using  an  image  plate  is  shown 
in  Fig.  2b  below  the  diagram  of  the  spectrometer.  The 
energy  range  is  8-12  keV,  encompassing  tungsten’s  La, 
Lp  and  Ly  lines.  The  LP2  line  (Fig.  3b)  is  shifted  to 
higher  energy  by  21  eV,  and  broadened  relative  to  the 
other  L  lines.  The  energy  shift  is  attributed  to  an  average 
ionization  of  +17  and  the  LP2  line  shape  is  attributed  to 
ionization  states  in  the  +10  to  +28  range  that  are  present  at 
different  times  during  the  pulse. 

C.  X-ray  distribution  measurements 

The  rolled  edge  technique  (Fig.  2c)  is  a  standard  way 
to  measure  radiographic  source  distributions,  usually  end- 
on  and  time-integrated.  The  PFRP  x-ray  source 
distributions  are  measured  side-on  with  -  1  ns  time 
resolution  using  an  array  of  20  pin  diodes.  The  array  is 
aligned  either  in  the  axial  (parallel  to  the  axis)  or  vertical 
(perpendicular  to  the  axis)  direction.  Each  1x1  mm  pin 
diode  is  part  of  a  linear  strip  of  10  diodes.  Two  parallel 
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Figure  3.  (a)  Temperature  vs.  time  (b)  LP2  line  shape  (squares,  time  integrated)  is  shifted  to  20  eV  higher  energy  and 
broadened  by  multiple  ionization  states  between  +10  and  +28  (from  Ref.  6). 


strips  are  staggered  to  make  an  array  of  20  diodes  with  0.5 
mm  between  adjacent  centers.  The  pin  diodes  are  cross- 
calibrated  and  time-synchronized  using  data  from  a  flat- 
field  shot  without  the  rolled  edge.  The  magnification  M  = 
(rolled  edge  -  pin  array)/(rolled  edge  -  PFRP)  was  2  for 
the  axial  measurements  and  3  for  the  vertical 
measurements. 

At  a  given  time,  the  signal  amplitudes  are  plotted  vs. 
position  (the  “edge  spread”  function).  For  a  Gaussian 


source  distribution,  -  exp-(r/a)^,  the  edge  spread  is  an 
offset  error  ftinction,  -  l+erf[(y-yo)/a].  Error-ftinction  fits 
to  the  edge  spreads  are  excellent  for  the  vertical  array  data 
and  good  for  the  axial  array  data.  The  fits  give  the  1/e 
widths,  ay  for  the  vertical  array  data  and  az  for  the  axial 
data. 

The  resulting  vertical  and  axial,  Gaussian  x-ray 
distributions  are  plotted  in  Fig.  4  at  several  times.  Above 
each  plot  is  a  visible-light  photograph  of  the  rod  taken 
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Figure  4.  X-ray  distributions  from  (a)  vertical  pin  diode  array,  (b)  axial  pin  diode  array. 
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Time  (ns)  after  start  of  x-rays 

Figure  5.  Mass  density  on  axis  vs.  time,  calculated  assuming  two  different 
rod  shapes. 


prior  to  the  shot.  A  blowup  of  the  rod  tip  is  superimposed 
on  each  plot  to  compare  with  the  x-ray  distributions.  At  t 
=  1  ns,  the  vertical  x-ray  distribution  (Fig.  4a)  extends 
well  beyond  the  original  rod  boundary  implying  the 
tungsten  is  expanded  prior  to  the  start  of  the  x-rays, 
possibly  the  result  of  joule  heating  by  the  generator 
current  flowing  through  the  small-radius  tip  prior  to  the 


Figure  6.  (a)  Vertical  scale  length  and  (b)  mass  per  unit 
length  as  hinctions  on  time.  Filled  symbols  are  at  t  =  1, 
2,  5,  10,  and  20  ns,  the  times  plotted  in  Fig.  4. 


formation  of  the  electron  beam.  From 
1-5  ns,  the  vertical  distribution 
becomes  more  narrow  (ay  decreases). 
The  axial  distribution  (Fig.  4b) 
originates  near  the  rod  tip,  increasing  in 
amplitude  from  1-5  ns.  From  5-20  ns 
the  axial  and  vertical  scale  lengths  both 
increase. 


III.  DATA  ANALYSIS 


The  tungsten  radial  mass  distribution 
is  assumed  to  be  the  same  as  the 
measured  vertical  x-ray  distribution. 
This  assumption  is  valid  if  the  electron 
energies  are  high  enough  so  their  range 
is  much  greater  than  the  plasma 
thickness,  and  the  plasma  is  transparent 
to  the  x-rays  recorded  by  the  pin 
diodes.  Both  of  these  conditions  should 
be  true  for  the  PFRP  experiments  on  Gamble  II.  The 
electron  range  exceeds  the  -  .05  g/cm^  tungsten  thickness 
for  energies  greater  than  0.15  MeV.  The  pin  diodes 
respond  to  x-rays  that  are  transmitted  through  a  13  mm 
thick  plastic  vacuum  window;  half  of  the  pin  diode  signal 
is  from  x-ray  energies  greater  than  300  keV.  The  tungsten 
HEDP  is  highly  transparent  to  these  x-rays.  The  axial  x- 
ray  distribution  and  the  a  priori  rod  photograph  are  used 
to  infer  the  HEDP  mass.  Together,  these  measurements 
are  used  to  calculate  the  radial  mass  distribution  and  the 
plasma  pressure. 

A.  Mass  density 

The  length  of  the  HEDP  is  taken  to  be  1.9  times  the 
axial  scale  length,  az.  This  choice  includes  82%  of  the 
radiation  and  discounts  the  lower-radiation  wings  of  the 
distribution.  The  rod  photographs  (Fig.  4)  are  used  to 
calculate  the  initial  tungsten  mass  between  the  tip  and  z  = 
1.9az.  The  mass  per  unit  length,  p,  is  assumed  to  be  this 
mass  divided  by  1.9az  for  the  tungsten  HEDP. 

The  radial  mass-density  distribution,  p(r),  is  derived 
from  the  vertical  x-ray  distribution  and  p.  Since  the 
vertical  distribution  is  Gaussian,  the  radial  distribution  is 
also  Gaussian  with  the  same  scale  length  ay.  The  mass 
distribution  is  then: 

p(r)  =  poexp-(r/ay , 
where  po  =  p/(7iay^). 

The  density  on  axis,  po,  for  each  rod  in  Fig.  4  is 
plotted  vs.  time  in  Fig.  5.  The  rod  shapes  were  not 
identical;  the  end  diameters  were  0.55  mm  for  shot  10240 
with  the  axial  array  and  0.75  mm  for  shot  10238  with  the 
vertical  array.  The  peak  densities  are  0.7  and  0.95  g/cm^ 
at  t  =  7  ns.  (This  uncertainty  would  be  eliminated  by 
measuring  both  x-ray  distributions  on  the  same  shot.)  At 
the  time  of  maximum  pressure  (4.5  ns,  see  below)  the 
density  is  0.45-0.75  g/cm^. 


319 


Figure  7.  Pressure  calculated  from  mass  dynamics.  Total  pressure  (black 
line)  is  sum  of  three  terms,  mass-acceleration  (red),  m-dot-v  (blue)  and 
magnetic  pressure  (green). 

The  ion  coulomb  coupling  parameter,  F,  is  the  ratio 
of  potential  energy  to  kinetic  energy  [7] : 


dominated  by  the  mass-time-acceleration 
term.  At  t  =  4.5  ns  there  is  a  local 
maximum  in  the  pressure  of  16  Mb.  The 
pressure  is  calculated  using  the  mass  per 
unit  length  for  the  rod  with  the  larger  tip 
diameter.  The  pressure  is  about  9  Mb 
when  the  other  rod  tip  is  used  in  the 
calculation.  This  large  uncertainty  could 
be  eliminated  by  measuring  both  x-ray 
distributions  during  the  same  shot. 

This  calculated  pressure  is  compared 
with  UekT,  using  the  mass  density  at  t  =  4.5 
ns  in  Fig.  5.  The  density  range  at  this  time 
is  0.45-0.75  g/cm^  or  n;  =  1. 5-2.5  x  10^' 
cm'^.  Assuming  (1+Z)  =  20,  He  =  (1+Z)  nj 
=  3-5  X  10^^  cm'^  and  nekT  =  1. 4-2.4  Mb. 
This  pressure  calculated  from  the 
expansion  of  the  plasma  is  higher  than 
these  values  by  about  a  factor  of  7.  The 
reason  for  the  higher  pressure  is  the  strong 
coupling  between  plasma  ions,  and  the 
inability  of  the  plasma  electrons  to 
completely  screen  the  ion-ion  coulomb 
forces.  The  number  of  electrons  in  a 
Debye  sphere  is  approximately  1.6,  so 
normal  screening  approximations  are  not 
valid. 
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Using  the  measured  PFRP  parameters:  Z  =  17,  A  =  184, 
pm  =  0.7,  T  =  30  eV,  gives  F  =  35,  therefore  the  PFRP 
plasma  is  strongly- coupled  (F  >  1). 


B.  Plasma  pressure 

The  expansion  of  the  HEDP  mass  distribution  is 
attributed  to  its  internal  pressure.  The  pressure  can  be 
calculated  from  Newton’s  second  law,  written  in  terms  of 
the  measured  parameters  of  the  x-ray  distributions: 

d^a  dudaA 

A - T^  +  —  — ^ 

dr  dt  dt  200m  ^ 

V  J  y 

The  first  term  in  brackets  is  the  mass-times-acceleration 
term;  the  second  term  is  the  mass-derivative  times 
velocity.  The  last  term  is  the  maximum  magnetic 
pressure.  The  time-dependent  parameters  ay  and  p  are 
plotted  in  Fig.  6a  and  b.  The  measured  values  (symbols) 
are  fit  by  polynomial  functions  that  are  used  to  calculate 
the  time-derivatives  in  the  pressure  equation. 

An  example  pressure  calculation  is  shown  in  Fig.  7. 
The  three  components  in  the  pressure  equation  are  plotted, 
with  the  total  pressure  indicated  by  the  black  line.  The 
initial  pressure  is  high,  but  this  is  probably  an  artifact  of 
electron  beam  motion  instead  of  mass  motion  (dp/dt  and 
day/dt  are  both  negative).  For  t  =  2-10  ns,  the  pressure  is 


2m, 


IV.  SUMMARY 

This  paper  reports  measurements  of  HEDP  parameters 
created  at  the  tip  of  a  tapered  tungsten  rod  in  a  PFRP 
diode.  The  plasma  parameters  are  in  the  Warm  Dense 
Matter  (WDM)  regime:  T  =  30  eV,  p  =  0.7  g/cm^  Z  =  17, 
and  P  =  16  Mb.  The  plasma  ions  are  strongly  coupled, 
with  F  =  35.  The  pressure  is  calculated  from  the 
measured  plasma  expansion,  independent  of  temperature 
and  ionization,  and  is  found  to  exceed  nkT  by  about  a 
factor  of  7.  This  HEDP  is  produced  by  a  “standard” 
pulsed  power  generator,  representing  a  new  way  to 
produce  WDM  in  the  laboratory. 

The  HEDP  parameters  presented  here  are  the  result  of 
measurements  made  during  different  shots.  All  the 
diagnostics  could  be  implemented  on  the  same  shot  in 
future  experiments  to  eliminate  uncertainties  from  shot- 
to-shot  variations  (due  to  different  rod  dimensions,  for 
example).  The  spectroscopic  diagnostic  could  be  made 
time-resolved  using  an  x-ray  streak  or  framing  camera. 
The  x-ray  distribution  measurements  could  be  improved, 
with  precise  alignment  between  pin  diodes  and  the  rod  tip, 
and  with  two  arrays  and  two  rolled  edges.  Two- 
dimensional  x-ray  imaging  [8]  would  be  useful  for 
comparison  with  2D  models.  An  improved  diagnostic 
suite  could  then  be  used  to  explore  the  dependence  of  the 
HEDP  parameters  on  variations  in  the  diameter  of  the  rod 
tip  or  the  pulsed  power  parameters. 
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